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O , Abstract 

| 7 j ; We find that in "two-photon" -like processes in the scalar (p 3 E model and also in 

\q • hadron-pair production arising from the collisions of a real (transversely polarized) 

and a highly virtual, longitudinally polarized, photon in QCD, there is duality be- 
tween two distinct nonperturbative mechanisms. These two mechanisms, one involv- 
ing a twist-3 Generalized Distribution Amplitude, the other employing a leading- 
twist Transition Distribution Amplitude, are associated with different regimes of 
factorization. In the kinematical region, where the two mechanisms overlap, duality 
is observed for the scalar tp\ model, while in the QCD case the appearance of duality 
turns out to be sensitive to the particular nonperturbative model applied and can, 
therefore, be used as a tool for selecting the most appropriate one. 
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1 Introduction 



The only known method today of applying QCD in a rigorous way is based on the factor- 
ization of the dynamics and the isolation of a short-distance part that becomes accessible to 
perturbative techniques of quantum field theory (see, [H [2j [3] and for a review, for instance, 



^ ■ [1] and references cited therein). Then, the conventional systematic way of dealing with 
the long-distance part is to parametrize it in terms of matrix elements of quark and gluon 
operators between hadronic states (or the vacuum). These matrix elements stem from non- 
perturbative phenomena and have to be either extracted from experiment or be determined 
on the lattice. In many phenomenological applications they are usually modeled in terms of 
various nonperturbative methods or models. 

Generically, the application of QCD to hadronic processes involves the consideration of 
hard parton subprocesses and (unknown) nonperturbative functions to describe binding ef- 
fects. Prominent examples are hard exclusive hadronic processes which involve hadron dis- 
tribution amplitudes (DAs), generalized distribution amplitudes (GDAs), and generalized 
parton distributions (GPDs) [21 [6], [TJ [8] . Applying such a framework, collisions of a real and 
a highly- virtual photon provide a useful tool for studying a variety of fundamental aspects of 
QCD. 

Recently, nonperturbative quantities of a new kind were introduced — transition distri- 
bution amplitudes (TDAs) [9l [101 [11] — which are closely related to the GPDs. In contrast 
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to the GDAs, the TDAs appear in the factorization procedure when the Mandelstam vari- 
able s is of the same order of magnitude as the large photon virtuality Q 2 , while t is rather 
small. Remarkably, there exists a reaction where both amplitude types, GDAs and TDAs, 
can overlap. This can happen in the fusion of a real and transversely polarized photon with 
a highly- virtual longitudinally polarized photon, giving rise to a final state which comprises 
a pair of pions. The key feature of this reaction is that it can potentially follow either path: 
proceed via twist-3 GDAs, or go through the leading-twist TDAs, as illustrated in Fig. CD 
Such an antagonism of alternative factorization mechanisms in this reaction seems extremely 
interesting both theoretically and phenomenologically and deserves to be studied in detail. 

The intimate relation between these two mechanisms in the production of a vector-meson 
pair was analyzed in [12] and it was found that these mechanisms can be selected by means 
of the different polarizations of the initial-state photon. In contrast, for (pseudo) scalar par- 
ticles, such as the pions, this effect is absent enabling us to access the overlap region of both 
mechanisms and their duality as opposed to their additivity. 

In this talk, we will report on the possibility for duality between these antagonistic mech- 
anisms of factorization, associated either with GDAs or with TDAs, in the regime where both 
Mandelstam variables s and t are rather small compared to the large photon virtuality Q 2 . 




Figure 1: Two ways of factorization: via the GDA mechanism and via the TDA mechanism. 



2 Regimes of Factorization within the ^-model 

Consider first the factorization of the scalar Lp\ model in Euclidean space. To study the 
four-particle amplitude in detail, it is particularly useful to employ the a-representation — see 
[7j. Then, the contribution of the leading "box" diagram can be written as (while details can 
be found in [T3] ) 
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where m 2 serves as a infrared (IR) regulator, s > 0, t > are the Mandelstam variables in 
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the Euclidean region, and D = a i- Assuming that q 2 = Q 2 is large compared to the mass 

i=l 

scale m 2 (which simulates here the typical scale of soft interactions) , the amplitude can 
indeed be factorized. As regards the other two kinematic variables s and t, one can identify 
three distinct regimes of factorization: (a) s < Q 2 while t is of order Q 2 ; (b) t <ti Q 2 while s 
is of order Q 2 ; (c) s, t <£i Q 2 . 
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Figure 2: The ratios R\ and R 2 as functions of s/Q 2 . 



Regime (a): The process is going through the s-channel. In this regime, the main 
contribution in the integral in Eq. (pQ) arises from the integration over a x when ot\ ~ 0: 
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Schematically this means that the propagator, parametrized by cti, can be associated with 
the partonic (hard) subprocesses, while the remaining propagator constitutes the soft part of 
the considered amplitude, i.e., the scalar version of the GDA. 

Regime (b): Here we have to eliminate from the exponential in Eq. ([T]) the variables 
Q 2 and s, which are large. This can be achieved by integrating over the region a 2 ~ 0. 
Performing similar manipulations as in regime (a), we find that the scalar TDA amplitude 
can be related to the scalar GDA via ^4tda( s 5 m2 ) = ^gda(^ s > rn 2 ). 

Regime (c): The relevant regime to investigate duality is when it happens that both 
variables s and t are simultaneously small compared to Q 2 , i.e., when s, t <C Q 2 . In this 
case, there are two possibilities to extract the leading Q 2 -asymptotics, notably, we can either 
integrate over the region a± ~ 0, or integrate instead over the region a 2 ~ 0. Clearly, these 
two options can be associated with (i) the GDA mechanism of factorization with the meson 
pair scattered at a small angle in its center-of-mass system or, alternatively, (ii) with the TDA 
mechanism of factorization. We stress that we may face double counting when naively adding 
these two contributions. We interpret such a behavior as a signal of an ingrained tendency 
for duality between the GDA (s-channel) and the TDA (t-channel) factorization mechanisms. 

In order to verify the appearance of duality we carry out a numerical investigation of 
the exact and the asymptotic amplitudes. In doing so, we introduce the following ratios 
R\ = A^ DA /A and R 2 = -A.q DA /A. Appealing to the symmetry of these ratios under the 
exchange of the variables s t, we take t/Q 2 to be 0.01 and look for the variation of the 
ratios with s/Q 2 . This variation is illustrated in Fig. [2] from which one sees that in the region 
where s/Q 2 is rather small, i.e., in the range (0.01, 0.05), both asymptotic formulae are 
describing the exact amplitude with an accuracy of more than 90%. This behavior supports 
the conclusion that, when both Mandelstam variables s/Q 2 and t/Q 2 assume values in the 
wide interval (0.001, 0.7), duality between the TDA and the GDA factorization mechanisms 
emerges. 



3 TDA- and GDA-Factorizations for 77* — > tttt 

Having discussed the appearance of duality between the GDA and the TDA factorization 
schemes within a toy model, we now turn attention to real QCD. To analyze duality, we 
consider the exclusive tt + 7t~ production in a 7t7l collision, where the virtual photon with a 
large virtuality Q 2 is longitudinally polarized, whereas the other one is quasi real and trans- 
versely polarized. Notice that the GDA and the TDA regimes correspond to the same helicity 
amplitudes. Given that the considered process involves a longitudinally and a transversally 
polarized photon, we are actually dealing with twist-3 GDAs |14J . On the other hand, for 
the twist-2 contribution, related to the meson DA, we use the standard parametrization of 
the 7r + -to- vacuum matrix element which involves a bilocal axial- vector quark operator pQ. 
Finally, the 7 — > 7r~ axial- vector matrix elements can be parametrized in the form, cf. [10J, 

(n-(pM(-z/2) lal5 [-z/2;z/2}^(z/2)\ 1 ( q ',e')) £ -f e' T ■ A T P a Ai(ar, £, t) , (3) 
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where P = (p2 + Q r )/2, and A = p 2 — q', and noticing that the symbol = means Fourier 
transformation and that the vector matrix element does not contribute here. To normalize 
the axial- vector TDA, A%, we express it in terms of the axial- vector form factor measured 
in the weak decay 7r — > luij [T5J dSl H3]. The helicity amplitude associated with the TDA 
mechanism reads = jF TDA e'^) • A T /|<f | with 

- M<|^|w,W,.) - ^> (4) 

where we have employed the 1-loop a s (Q 2 ) in the MS-scheme with Aqcd = 0.312 GeV for 
Nf = 3 [17]. [Note that there is only a mild dependence on Aq CD .] 

Turning now to the helicity amplitude, which includes the twist-3 GDA, we anticipate 
that it can be written as (see, for example, [H]) A?^J\ = J^^e'O) . A T /\q \ with 
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where the partial derivative is defined by = d/d(2( — 1). In deriving ([5]), we have used for 
the twist-3 contribution the Wandzura-Wilczek approximation. Duality between expressions 
(jl]) and (jSJ) may occur in that regime, where both s and t are simultaneously much smaller 
in comparison to the large photon virtuality Q 2 . More insight into the relative weight of the 
amplitudes with TDA or GDA contributions can be gained once we have modeled these non- 
perturbative quantities. We commence with the TDAs and, assuming a factorizing ansatz 
for the t-dependence of the TDAs, we write Ai(x, £, t) = 2-^- F^(t)Ax(x, £), where the t- 
independent function Ai(x,£) is normalized to unity. To satisfy the unity-normalization 

condition, we introduce a TDA defined by A x {x,l) = Af n - norm (x, 1)/ ) dxA™ n ~ m (x, 1) 

-1 

and continue with the discussion of the t-independent TDAs. Recalling that we are mainly 
interested in TDAs in the region £ = 1 [1, 2J, it is useful to adopt the following parametrization 

A non-norm (X; X ) = (1 _ ^ ^ + axC ^ 2 \ X ) + a 2 Cf /2] \x) + a A Cf /2 \x)^ , (6) 

where a±, a 2 , a 4 are free adjustable parameters, encoding nonperturbative input, and the 
standard notations for Gegenbauer polynomials are used. It is not difficult to show that the 




Figure 3: Helicity amplitudes jF TDA and jF GDA as functions of Q 2 , using a\ = —0.5 found in 
lattice simulations. The value of s/Q 2 varies in the interval [0.06, 0.3]. 

TDA expressed by Eq. ([6]) results from summing a D-term, i.e., the term with the coefficient 
ax, and meson- DA-like contributions. For our analysis, we suppose that a% = do [8], which is 
equal to —0.5 in lattice simulations. With respect to the parameters a 2 and 04, we allow them 
to vary in quite broad intervals, notably, a 2 G [0.3, 0.6] and 04 G [0.4, 0.8], that would cover 
vector-meson DAs with very different profiles at a normalization scale /i 2 ~ 1 GeV 2 (see, for 
example, [IE])- The function C) is rather standard and well-known (details in [5| [T3]). 

We close this section by summarizing our numerical analysis of [IB]. We calculated both 
functions jF TDA and JF GDA , and show the results in Fig. [3j The dashed line corresponds to 
the function JF TDA , where we have adjusted the free parameters to a 2 = 0.6, a 4 = 0.8. The 
results, obtained for rather small values of these parameters, are displayed by the broken lines 
in the same figure. The dotted line denotes the function jF TDA with a 2 = 0.5 and 04 = 0.6, 
whereas the dash-dotted line employs a 2 = 0.3 and a 4 = 0.4. For comparison, we also include 
the results for JF GDA . In that latter case, the dense-dotted line corresponds to the GDA 
amplitude, where the expression for B\ 2 has been estimated via Eq. (20) of [13], while the 
solid line represents the simplest ansatz for B\ 2 with R n = 0.5. From this figure one may 
infer that when the parameter i? 12 , which parametrizes the GDA contribution, is estimated 
with the aid of the Breit-Wigner formula (provided s, t <C Q 2 ), there is duality between the 
GDA and the TDA factorization mechanisms. Hence, the model for $1(2, (), which takes into 
account the corresponding resonances, can be selected by duality. 



4 Conclusions 

We have provided evidence that when both Mandelstam variables s and t turn out to be 
much less than the large momentum scale Q 2 , with the variables s/Q 2 and t/Q 2 varying in 
the interval [0.001, 0.7], then the TDA and the GDA factorization mechanisms are equivalent 
to each other and operate in parallel. We have also demonstrated that duality may serve as 
a tool for selecting suitable models for the nonperturbative ingredients of various exclusive 
amplitudes entering QCD factorization. In this context, we observed that twist-3 GDAs 
appear to be dual to the convolutions of leading-twist TDAs and DAs, multiplied by a QCD 
effective coupling. 
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